Practionalization in Josephson Junction Arrays Hinged by Quantum Spin Hall edges 
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We study a novel superconductor-ferromagnet-superconductor (SC-FM-SC) Josephson array de- 
posited on top of a two-dimensional quantum spin Hall insulator. The Majorana bound state at 
the interface between SC and FM leads to charge-e tunnelling between neighboring superconductor 
islands, in addition to the usual charge-2e Cooper pair tunnelling. Moreover, because Majorana 
fermions encode the information of charge number parity, an exact Z2 gauge structure naturally 
emerges and leads to many new phases, including a deconfined phase where electrons fractionalize 
into charge-e bosons and topological defects. A new SC-insulator transition has also been found. 
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I. INTRODUCTION 

Superconductor Josephson junction arrays have been 
studied extensively in the past. At low temperature, 
the phase of a superconductor island becomes a quan- 
tum degree of freedom, which is conjugate to the Cooper 
pair density. Competition between charging energy and 
Josephson tunnelling leads to a superconductor-insulator 
transition at zero temperature, which is usually stud- 
ied using quantum rotor models. Such a boson-only ap- 
proach is no longer adequate if low-energy quasi-particles 
are present. Recently, it was proposed that if a s-wave 
superconductor (SC) and ferromagnet (FM) junction is 
hinged by the edge of a quantum spin Hall (QSH) insula- 
tor discovered recently^, a zero-energy Majorana bound 
state will be localized at the SC-FM interface 2 . Al- 
though the Majorana fermions do not carry the global 
U(l) charge of electrons, they do encode the informa- 
tion of fermion parity i.e. the even-odd of the electron 
number, hence Majorana fermion will participate in the 
tunnelling of charges between two SC islands, and enrich 
the physics of the Josephson array. 

So far most work on topological insulators have been 
focusing on the band structure or weak interaction ef- 
fects. Recently strong correlation effects for topological 
insulators have attracted more and more attentions^—. 
Although it has been proposed that quantum computa- 
tion may be realized based on Majorana bound states 
localized in topological defects 3 , very little was studied 
about the many-body or strong correlation effects that 
the Majorana fermions can participate in. In this work 
we will focus on the correlation physics in the Majorana 
fermion assisted Josephson array, and we will show that 
such Josephson arrays have unusual fractionalized phases 
and phase transitions. 



II. ONE DIMENSIONAL CASE 

Let us warm up with a simple one dimensional ge- 
ometry, as depicted in Fig. QJi. We denote the location 
of every SC island by coordinate j, and denote the two 



Majorana fermions around each SC island as 7^1 and 
7^2- The FM islands have uniform magnetizations that 
are perpendicular to the angular momentum carried by 
the electrons of the QSH edge states, which opens up 
a gap for the QSH edge states. As was pointed out by 
Refii, these two Majorana fermions correlate with the 
fermion number on this SC island through the constraint 
*7j,i7j.2 = (— 1)™j. The tunnelling between SC islands 
has to be consistent with this constraint, and the follow- 
ing two terms are allowed 8 : 

H a = }^ -h cos(<f>j - <t>j+ x ), 
3 

H t 2 = ^ -tiiljalj+x.i cos[(^j - (j) j+x )/2\. (1) 
3 

Here <f>j is the phase angle of the SC island j, and e^ j 
increases the electron number rij by 2. Hti is the ordi- 
nary Josephson tunnelling term, and H t 2 is the charge-e 
tunnelling assisted by the Majorana fermions, which is 
now allowed because with the Majorana fermions there 
is no longer a Cooper pair breaking gap between even 
and odd electron number on each SC island^ i.e. electron 
can reside across the SC island nonlocally through ma- 
jorana zero modes 7^1 and 7^,2 • Since H t \ is a second 
order effect that involves a Cooper pair breaking inter- 
mediate state, in the limit with dominant Cooper pair 
energy, H t \ is ignorable. Inclusion of small H t \ will only 
quantitatively change the physics discussed in this paper. 

It is well-known that the one dimensional Majorana 
fermion is equivalent to a transverse field quantum Ising 
model, and the Ising variables are defined on the links of 
the Id lattice denoted as + x) in Fig. [TJ 
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jj+x = II ir YkA7k,2, cr z ] j+x = i-Yj,2lj+xA- (2) 
k<j 



Now the full Hamiltonian can be written as 

H = u ( n i ~ «) 2 " t^h+x cos (y - - 
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which is subject to the constraint 

"j-x,j a j,3+x{ L ) — ' 



fj+x ■ 
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FIG. 1: (a), the lattice structure for case 1. The green circles 
are the Majorana fermions, and the dashed brown line is the 
edge of the underlying QSH insulator, (b), the phase diagram 
for ordinary SC Josephson array, (c), the phase diagram for 
Josephson array in our case, where the number of MI lobes 
doubles compared with (b). 



The U term represents a charging energy. The Hamilto- 
nian Eq. [3] takes the standard form of a Z 2 gauge field a z 
coupled to matter field </>/2. For this ID system the Z-i 
gauge field cr? a +x can be absorbed into the rotor variable 
<j)j through the following duality mapping: 



2 



1 - 7? 



rf.(5) 



Tj,Tj = ±1 are Ising operators defined on SC islands, 
and they satisfy the algebra of Pauli matrices. 8j and nj 
satisfy the rotor phase-number algebra: [9j,nk] = iSjk- 
Now this model can be written with the new variables as 
ordinary Bose-Hubbard model^ in ID: 



H = 



3 



U{iij — n) 2 — t 2 cos( 



(6) 



and the gauge constraint operator ^jj-x^j+xi^^Jj = 
tJ{— l) nj commutes with e l0j . 

The phase diagram of the model Eq. [3] is identical to 
model Eq. |5J with integer n, there is a Mott-insulator 
(MI) phase with fixed fermion number on every SC is- 
land when tijU <C 1, and a superfluid (SF) phase when 
tijU ^> 1. The phase diagram is shown in Fig.[T]:. When 
n is precisely half-integer, since in this model there is no 
nearest neighbor density repulsion, the MI phase van- 
ishes. Compared with the ordinary Josephson array, the 
number of MI lobes is doubled in our case (Fig.QJ, c), ba- 
sically because there is no charge gap between even and 
odd particle filling. The MI phase is nondegenerate, ev- 
ery SC island has fixed charge number. In the SF phase, 
besides the Cooper pair operator e 1 ^, the following string 
operator also has algebraic correlation: 
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-i(j> j+L /2 



>~(6 



^>~(1/L)*.(7) 



with the Luttinger parameter K tuned by tijU . 

Single electron can be injected into the system through 
the Majorana fermion bound state. One legitimate rep- 
resentation of electron operator is Cj\ ~ e I ^/ 2 7j j i— 
The spin index does not appear in the electron operator 
because the spin degeneracy is lifted by FM islands and 
the spin-orbital couplings in underlying QSH edge states. 
Also, under transformation <j)j — > <pj + 27r, the solution 
of the Majorana bound state -fj, a changes sign, there- 
fore the physical electron operator remains invariant. In 
terms of the bosonic variables, the electron operator can 
be expressed as 



n 

k<j 



(J k-x,k a j-x,j e 



exp[i7r n,k]e l 
k<j 



(8) 



The a x in the product Eq. [S] guarantees the fermionic 
statistics between two electron operators. The correla- 
tion function between two electron operators is 



L-l 

exp[l7T '/;/, f 

k=j-L 
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The Bosonic representation of electron operator Eq. [8] 
takes exactly the same form as the standard fermioniza- 
tion of the Bose variables 9 and n in one dimension. 

The phase transition at the integer filling is a 
Kosterlitz-Thouless (KT) transition at K — 1/4, which 
physically corresponds to proliferating 27r kinks of 9 de- 
fined in Eq. which is equivalent to a 47r kink of <j). 
At this transition, the scaling dimension of Cooper pair 
operator is 1/2, while in ordinary ID Josephson array, at 
the KT transition the Cooper pair operator has scaling 
dimension 1/8, because normally the transition is driven 
by the proliferation of 2u kink of 0. The difference be- 
tween these two cases stems from the charge-e tunnelling 
enabled by the Majorana fermions. In terms of 9 intro- 
duced in Eq. [SJ H t \ is simply Yl j_ — ti cos(29j ~ 29j +x ), 
therefore turning on H t \ in Eq. |6] will not change the 
phase diagram, as long as there is no pairing gap be- 
tween even and odd filling. The random one-dimensional 
lattice of Majorana fermions at the edge of a topologi- 
cal insulator, which is similarly connected to the random 
Ising model, has been discusse d 10 ' 11 . 



III. TWO DIMENSION, REGULAR 
STRUCTURE 



Now let us move on to the 2D structure as depicted in 
Fig.[2Ji- We place the SC and FM islands on top of a QSH 
insulator with commensurate holes. The internal edges 
of these holes hinge all the islands in this lattice. We first 
assume that the FM island is very thin compared with SC 
island, therefore the intra-island tunnelling between Ma- 
jorana fermions is negligible compared with inter-island 
tunnellings. We denote every SC island as site j, and 
denote the four Majorana fermions around each SC is- 
land as (j, a) with a — 1 • • ■ 4. The effective lattice is 
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shown in Fig. [2}). Now the gauge constraint becomes 
7j,i7i,27j.37j,4(^l)' li = 1- Again we can map the Majo- 
rana fermions to Zi gauge field as following: 

a lj+v = II *7fc,47fc,3, (Tj t j +X = ]J *7fc,i7fc,2- (10) 

The entire Hamiltonian can be parallelly generalized from 
its Id counterpart Eq. [3] 

H = Y J U{n 3 - nf - £ t 2 al j+v cos(^ - ^)(11) 
Again this Hamiltonian is subject to the gauge constraint 

-X x x x ( _~\\ n j 

a j,j-x"j,j+x cr j,j-y a j,j+y - \ 1 ) ■ l lz 7 

The 2D Z 2 gauge field is drastically different from ID, 
in the sense that it has a nontrivial liquid phase even 
when (f>j is disordered. In the MI phase of SC islands, in- 
tegrating out the gapped fluctuation <pj will induce gauge 
invariant dynamics for Zi gauge field: 

Hring = y ^ ^-^- (7 j,j+x' 7 j+x,j+x+y a 'j+y,j+x+y' J j,j+yJ (13) 
3 

with K ~ t^/U 3 . This is a standard ring exchange term 
of Z 2 gauge field. This term favors the ring product of 
a z around each unit plaquette to be 1. In ordinary Z 2 
gauge field, this ring exchange term will compete with 
the Z2 string tension term v ~haf i+v , and when K 
dominates h the system is in a Z2 liquid phase with 
topological order which cannot be described by local or- 
der parameter—. When h dominates K, the system is 
in a confined phase without topological order, where Zi 
charged matter is not just gapped, but also confined spa- 
tially by a linear potential. Physically these two phases 
can be understood by the behavior of "vison", which is 
a topological excitation with the product !!□ ctZ = — 1 
on one plaquette. The vison carries a global Zi charge, 
because one can only create/anihilate a pair of vison by 
operator crf i+l/ . In the liquid phase, the vison number is 
conserved mod 2, while in the confined phase the global 
Z2 symmetry is spontaneously broken. This effect is 
manifested in the dual description of the Z2 gauge field, 
which is formulated through the mapping: Y\o a? ~ = T ^ 1 
(Tja,,, = r-f t z . j denotes the dual lattice shown in 

JJ+V J-x j J 

Fig. [2f). Therefore the ordinary quantum Z2 gauge field 
is dual to a 2d transverse field quantum Ising model 
H = y-.„ -hT- z T-* - Kt?. When K » h, this Ising 
model is in the disordered phase, where the Z2 conserva- 
tion of t x (vison number) is preserved; when /i> K , the 
Z2 global symmetry of t z is spontaneously broken. 

In our case, operator crJ J+l/ is a nonlocal product of 
Majorana fermions, hence the string tension term g^a+v 
cannot exist in the Hamiltonian. Therefore in the dis- 
ordered phase of <j>j (MI of SC islands), the local vi- 
son number commutes with the Hamiltonian i.e. vison 




FIG. 2: (a), the Josephoson lattice structure of case 2. (6), 
the effective lattice of the Josephson array in (a). The SC 
islands are denoted by j, and the Z2 gauge field a x and a z 
are defined on links + v). The dual Ising variables r a 
are defined on the sites j of the dual lattice, represented by 
red dashed lines, (c), case 3 with compressed y direction, 
the intra-island tunnelling Ji7j,37j,4 is denoted by the pink 
dashed lines, (d), the phase diagram of case 3 with integer n, 
plotted against tz/U and J/U. Phase A, B and C represent 
the SF, the Z2 liquid, and the a x ordered phase respectively. 
The green line between phase A and B is a 3d XY transition, 
and the transition between B and C is a first order transition 
described by Eq. 1181 The direct transition between phase A 
and C may expand into a stable roton liquid phase. 

is completely static. Hence the Z2 gauge field is in its 
liquid phase. The Z2 liquid phase is deconfined i.e. an 
extra electron will carry an infinite string of a x due to 
the gauge constraint, but the energy cost is finite. For 
instance, the electron operator at site j, 1 (Fig. [2]j) can 
be written as 

CjA „ e^'\, x ~ J] a z k _ Xik a x _ XJ e^/ 2 , (14) 
k<j 

where the product includes all the a z on x— links to the 
left or below site j (Fig. |3]i). creates a pair of 

vison excitations, and since in the Z2 liquid state the 
vison excitation is gapped and conserved, the least energy 
consuming behavior of an injected electron is to form 
a bound state with two visons, and become a charge-e 
boson represented by following operator: 

&i~n <7 *-».* ex p(^/2). (15) 

k<J 

The SF phase can also be viewed as a condensate of bj. 
If an electron is injected into this system, it will fraction- 
alize into a mobile charge-e boson and a pair of static vi- 
sons. If a Cooper pair is injected into this system, it will 
fractionalize into two bosons, instead of two electrons. 
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FIG. 3: (a), the electron at site j,l (Fig. (2}i) can be rep- 
resented as exp(i<j>j /2) times a string of cr z , denoted by red 
circles; By binding two visons (blue circles), the fermion be- 
comes a boson, (b), the Wilson loop operator as a product 
of <7 Z (red circles), in the dual formalism is a product of r x 
on the dual lattice (blue circles). In the Z2 liquid phase, the 
perturbation of J term will only change the Wilson loop at 
the corner, represented by the green dashed square. 



In this liquid phase, the boson and the vison will have 
mutual semion statistics i.e. when a boson bj encircles 
a vison through a close loop, the system wave-function 
acquires a minus sign-^. 

The SF phase has vortex with hc/2e magnetic flux, 
which is bound with a vison. If we start from the SF 
phase, the Z2 liquid phase can be viewed as a conden- 
sate of the double vortex with hc/e magnetic flux of SF 
phase, since the MI liquid phase still has vison conserva- 
tion, and the product of quantum circulation of the vor- 
tex and the charge in the vortex condensate is a constant: 
QvQe = ha^. The Z2 gauge field has no gapless photon 
excitations, hence at the transition between SF phase 
and the Z2 liquid phase, the Z2 gauge field does not in- 
troduce further anomalous dimension to ip ~ exp(i<p/2). 
For instance, when n is integer this transition is a 3d 
XY transition described by order parameter ip. However, 
the physical Cooper pair operator ip 2 ~ exp(i0) gains a 
rather large anomalous dimension, which has been cal- 
culated by various methods^: 77^,2 ~ 1.47. Similar situ- 
ations were discussed at the transition between Z2 spin 
liquid and spiral order in frustrated magnet o 15 ' 16 . 

There is another equivalent way of describing the frac- 
tionalized Z2 liquid phase i.e. the disordered phase of <fi. 
Let us assume the filling on every SC island is even, and 
we can define Ising variables 

Mi = *7j,i7i,3 = *7j,27j,4, 



Mi = *7j,37i,2 = Hj.MjA: 
Ml = i'7i,i7i,2 = «7i,47i,3- 



(16) 



\x a satisfy the algebra of Pauli matrices. Notice that un- 
like a a introduced before, now the Ising variables fi a are 
defined on the sites of the lattice instead of the links. 
Now the ring exchange term H r i ns reads 



This is precisely the model introduce in Refill, as a 
concrete example with topological orders. This model 
is equivalent to the toric code modeli 2 -, which becomes 
manifest once we switch the definition of [i x and ji y in 
Eq .1161 for one of the two sublattices of the square lattice. 



IV. TWO DIMENSION WITH INTRA-ISLAND 
TUNNELING 

Now we consider a lattice structure slightly different 
from the previous section, with the entire system com- 
pressed in the y direction, until there is a considerable 
intra-island tunnelling — Ji"/j,3"fj,4- Using the mapping 
derived in last section, this operator is J a j,j- y a jj+ v 
in terms of the Z2 gauge variables, and this term be- 
comes a ring exchange on the dual lattice Ojj-yVjj+y = 
t*t-? t-. T-. . which represents a pair hopping of vi- 

J J— a: 3-x—y 3—y' r r rr a 

sons. Notice that t z are defined on the sites of the dual 
lattice instead of the links. The Hamiltonian of the dual 
Ising variables in the MI phase reads 



H 



- Jt-t- t- t- 

j 3—x j—x—y j-y 



(18) 



y 



(17) 



This is precisely the model studied in RefJ^ in the con- 
text of p ± ip superconductor Josephson array. The sym- 
metry of this model is quasi-local in the sense that we can 
reverse the sign of r z along any column or row arbitrarily, 
without changing the Hamiltonian. Physically this sym- 
metry means that the vison number has Z2 conservation 
along each row and column on the dual lattice. 

The phase diagram of model Eq. [18] has been studied 
both analytically and numerically^"^. Because of its 
special self-duality structure, it is expected that a tran- 
sition occurs at precisely J = K ~ t^/U 318 ^ 1 . When 
K > J, the system is in the Z2 liquid phase with topo- 
logical order; when J > K, the topological order van- 
ishes and the ground state is nondegenerate, although 
formally the nonlocal operator <j x has nonzero expecta- 
tion value. Due to the absence of string tension —hcr x j +u , 
Z2 charged matter is deconfined in both phases, which is 
very different from the ordinary Z2 gauge field. Mean 
field argument as well as numerical results suggest that 
the transition at J = K is first orde r 19 i 20 ' 22 . In the Z 2 
liquid phase, the J term enables the visons to move in 
pairs, therefore an injected electron will fractionalize into 
two mobile parts: charge-e boson and vison pair. Notice 
that unlike the ordinary Z 2 gauge field, this pair of visons 
cannot annihilate each other, due to the quasi local con- 
servation of vison numbers discussed in last paragraph. 

The difference between these two phases can be further 
characterized by the Wilson loop. In the phase with J > 
K, since the zeroth order ground state with K = is an 
eigenstate of t z , the Wilson loop (Ilc^) = Ilyt 7 " 2 can 
be calculated perturbatively with expansion of K , and it 
falls off according to an area law: (Y[ c <J Z ) ~ [K/J) — . 
Here C and A represent a closed loop and the area en- 
closed inside this loop respectively. In the Z2 liquid phase 
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with K > J, it is usually expected that with the presence 
of transverse field TV ~^ a j,j+u * ne Wilson loop falls off 
with a perimeter law, which can also be revealed with 
a perturbation of J/K on the ground state with J = 0, 
where the Wilson loop is a constant. However, In our 
situation with the intra-island tunnelling, the first order 
expansion of J term in Eq. [TH] will only change the Wil- 
son loop at the corners of loop C (Fig. Therefore we 
expect the Wilson loop in the Z 2 liquid phase falls off 
as a special "corner law" (Hc ^) ~ e~^ J l K , j\f is the 
number of corners of this Wilson loop C. 

Based on the analysis above, when J is small, by re- 
ducing tijV from infinity we will first drive a transition 
from the SF phase to a Zi liquid phase, and then enters 
an "area law" MI phase through a first order transition. 
When J is large enough, there can be a direct transition 
between the SF phase and the "area law" phase. This 
transition can be viewed as proliferating the hc/2e vor- 
tices of the SF phase which can only move in pairs due 
to the quasi-local conservation of visons. This type of 
paired directional vortex dynamics was the key of the 
roton liquid phase proposed befor o 21 ' 24 ' 25 , which is a sta- 
ble phase with gapless vison excitations and quasi one 
dimensional dispersions. Therefore the direct transition 
in Fig. might expand into a stable roton liquid phase. 
We will leave this possibility to future studies 2 ^. 

If we turn on not only —Jji'Jj,3'yj,4,, but also 
—Jji~fj,2jj,3 and — </^7j,37j,i, after introducing Ising 
variables /i a as Eq. [T6l the model describing the system 
becomes 

3 

••/;/'.; • -/f/'i • •/;/';• (19) 



It will be interesting to do a full analysis of all the 
possible phases of this model with different choices of 
site dependent transverse fields Jj 1 . For instance, with 
jx _ jy _ anc j jz j g s j^. e independent, this model 

reduces to the model in RefJ£. If J% = J V B ^ (A 
and B are two different sublattices of the square lattice) , 
while all the other transverse fields are zero, this model is 
equivalent to the toric code model with one component of 
transverse magnetic field, and there is a confine-deconfine 
phase transition driven by this transverse field. 



V. SUMMARY 

In summary, Josephson array with Majorana fermion 
zero modes around each SC island is described by a pre- 
cise Zi gauge field and matter field formalism, which 
leads to unusual fractionalization features in both one 
and two dimensions. Various pseudo-spin models with 
topological phases can be realized with the Josephson 
array. The fractionalization can be measured with single 
electron tunnelling experiments, because an electron will 
fractionalize into a boson and topological defects, and 
the single electron green's function becomes a convolu- 
tion of two fractional excitations. This will be discussed 
in more details in another paper. In our current paper 
we focus on the disordered phase of the Josephson array 
with fractional excitations, but when H t \ is nonzero or n 
is away from integer, many interesting phase transitions 
can occur inside the SF phase, we will also study these 
physics in future 2 ^. 
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